The effect of the three-spin interaction and the next-nearest neighbor interaction on 
the quenching dynamics of a transverse Ising model 
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We study the zero temperature quenching dynamics of various extensions of the transverse Ising 
model (TIM) when the transverse field is linearly quenched from — oo to +oo (or zero) at a finite 
and uniform rate. The rate of quenching is dictated by a characteristic scale given by r. The density 
of kinks produced in these extended models while crossing the quantum critical points during the 
quenching process is calculated using a many body generalization of the Landau-Zener transition 
theory. The density of kinks in the final state is found to decay as t'^^"^ . In the first model 
considered here, the transverse Ising Hamiltonian includes an additional ferromagnetic three spin 
interaction term of strength J3. For J3 < 0.5, the kink density is found to increase monotonically 
with J3 whereas it decreases with J3 for J3 > 0.5. The point J3 = 0.5 and the transverse field 
h — —0.5 is multicritical where the density shows a slower decay given by r~^^^. We also study 
the effect of ferromagnetic or antiferromagnetic next nearest neighbor (NNN) interactions on the 
dynamics of TIM under the same quenching scheme. In a mean field approximation, the transverse 
Ising Hamiltonians with NNN interactions are identical to the three spin Hamiltonian. The NNN 
interactions non-trivially modifies the dynamical behavior, for example an antiferromagnetic NNN 
interactions results to a larger number of kinks in the final state in comparison to the case when 
the NNN interaction is ferromagnetic. 

PACS numbers: 



I. INTRODUCTION 



The critical dynamics of classical systems have been 
studied extensively in last three decades while the study 
of the dynamics of a quantum system when swept across 
a quantum critical point (QCP) is fairly recent and not 
yet fully understood. The vanishing of energy gap be- 
tween the ground state and the first excited state of the 
quantum Hamiltonian signals the existence of a QCP— 1^. 
At a QCP, the correlation length as well as the relaxation 
time diverge, a phenomenon known as the critical slowing 
down. This diverging timescale makes it impossible for 
any system to cross the quantum critical point without 
excitations from the ground state. The dynamics there- 
fore is non-adiabatic in contrast to an adiabatic evolu- 
tion where the system sticks to the instantaneous ground 
state through out the quenching process. In recent years, 
there has been an upsurge in the study of dynamics close 
to a quantum critical point clearly indicating a growing 
interest in the field^^^^^i^^^^iii^i^. 

One of such attempts was to extend the Kibble's theory 
of defect production introduced to explain early universe 
behavior— to the second order quantum phase transi- 
tions. This method of calculating the density of defects 
is known as the Kibble-Zurek mechanism (KZM)^'^. The 
theory of KZM for a classical second order phase tran- 
sition is based on the universality of the critical slowing 
down and leads to the prediction that the linear dimen- 
sion of the ordered domains scales with the transition 
time T as where w is some combination of critical 
exponents. KZM has been confirmed by numerical simu- 
lations of time dependent Ginzburg-Landau modeU^ and 
also for various experimental system&iZ. On the other 



hand, for He-4 superfluid transition, KZM could not 
be verified experimentally^^. Hence more experiments 
are clearly needed to put KZM theory on a stronger 
footing. The same idea has been applied to study the 
dynamics across a zero temperature QCP by different 

The extended KZM for the zero-temperature quantum 
transitions relies on the fact that during the evolution 
when the system is close to the static critical point, the 
relaxation time diverges in a power-law fashion. The 
non-adiabatic effects become prominent when the time 
scale associated with the change of the Hamiltonian is of 
the order of the relaxation time. The loss of adiabatic- 
ity while crossing a quantum critical point can be quan- 
tified by estimating either the density of defects (e.g., 
the density of oppositely oriented spins in Ising mod- 
els) in the final state^'^'^'^'^ or the fidelity of the final 
state with respect to the ground statc^^ or the residual 
energ^ T^^i^^i^^i^^i^^i^"^ . The argument given above imme- 
diately leads to a (l/\/r)-dependence of the density of 
defects on the characteristic timescale r of the quench- 
ing. The residual energy is defined as the difference be- 
tween the energy of the evolved ground state and the 
true ground state. This residual energy for the inte- 
grable disorder free systems is trivially proportional to 
the density of kinks with the proportionality constant 
being equal to the strength of interaction. In an op- 
timization approach popularly known as the "quantum 
annealing" i9-20.gi, i22, 23^24 ^ ^j-^g strength of the quantum 

fluctuations is slowly reduced to zero so that a disor- 
dered and frustrated system of finite size is expected to 
reach adiabatically its true classical ground state. In 
the present literature, the expressions "annealing" and 
"quenching" are used synonymously. The residual en- 
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ergy turns out to be a more appropriate measure of non- 
adiabaticity for the annealing approach. In a recent work 
by Caneva et. al^, it has been shown that for a disor- 
dered quantum Ising spin chain, the residual energy and 
the density of kinks show different scaling behavior with 
T. Recently a general analysis has been carried out of 
the breakdown of the adiabatic limit in low-dimensional 
gapless systems?^. 

In this paper, we will concentrate on the estimation of 
density of defects produced during the dynamics of three 
different types of model Hamiltonians, all of them being 
exactly solved, at least in a mean field level, via the Jor- 
dan Wigner transformation. These three Hamiltonians 
are extensions of the TIM with an additional interac- 
tion term in each and our aim is to study the effect of 
such interactions on the density of defects produced dur- 
ing the quenching. The additional terms are i) a ferro- 
magnetic three spin interaction^^ ii) an antiferromagnetic 
next nearest neighbor interaction and iii) a ferromagnetic 
next nearest neighbor interaction, respectively. We con- 
sider the unitary evolution of the system prepared in the 
ground state of the initial Hamiltonian which crosses its 
equilibrium critical line as the system evolves. As de- 
scribed later, in all the cases, the fermionization of the 
Hamiltonian reduces it to a quadratic form and hence 
one can reduce the dynamics of a many-body Hamilto- 
nian effectively to a 2 x 2 Landau Zener problem'^* in the 
fourier representation. 

The paper is organized as follows. Section II includes 
a detailed discussion on the analytical diagonalization of 
the transverse Ising Hamiltonian with a three spin in- 
teraction term. In section III, we have described the 
transverse quenching scheme along with the results for 
the above model. We have presented a comparison be- 
tween the three spin Hamiltonian and the Hamiltonians 
with next nearest neighbor interactions when treated at 
a mean field level in section IV. A brief summary of the 
work is presented in the concluding section with a brief 
discussion based on the recent developments in this field. 



II. MODEL AND THE PHASE DIAGRAM 

The Hamiltonian of a one-dimensional three spin in- 
teracting transverse Ising system is given by21 

i 

+ J.5]«+J, (1) 

where and cr^ are non-commuting Pauli spin matrices, 
Jx is the strength of the nearest neighbor ferromagnetic 
interaction while J3 denotes the strength of the three 
spin interaction. In the limit J3 0, the model reduces 
to the celebrated transverse Ising model studied exten- 
sively in recent yearsii^. By a duality transformation"^ \ 
the above Hamiltonian can be mapped to a transverse 
XY model with competing (ferro-antiferro) interactions 



in the x and y components of the spinsi. Interestingly, 
even in the presence of the three spin interaction term, 
the Hamiltonian given by Eq. (1) is exactly solved by the 
Jordan- Wigner (JW) transformation32i^i^ which maps 
this interacting spin system to a system of noninteract- 
ing spinless fermions. Moreover, this three spin term is 
found to be irrelevant in determining the quantum crit- 
ical behavior of the system. The critical exponents are 
the same as that of Ising model in a transverse field ex- 
cept for the case Jx = ,and J3 = h. For the sake of 
completeness, let us now provide a brief discussion on the 
diagonalization of the Hamiltonian. 

In the JW-transformation, the Pauli matrices are 
transformed to a set of fermionic operators (c;) defined 
as 

i-l 

< = 24c, -1 (2) 

with cr^ = (cr"= + \a^)l2 and cr" = (cr^ - icr'')/2, and 
satisfy the standard anticommutation relations 

{c|, Cj\ = (5y , {c|, c]} = {q, Cj} = 0. 

We shall work in the basis in which is diagonal so that 
the presence of a fermion at a particular site i corresponds 
to an up spin (i.e., eigenvalue -1-1 of the operator erf) 
at that site. Using a periodic boundary condition, the 
Fourier transform of the Hamiltonian can be cast in the 
form 

U = ^[^(/^ + cos/c — J3Cos2/c)(cj,Cfe -I- c^j,c_fc) 

-t-i(sinfc - J3sin2A;)(4.cLj. + CkC-k)]- (3) 

Clearly, in the momentum representation of c-fermions, 
the Hamiltonian is quadratic and is translationally invari- 
ant. Using the Bogoliubov transformation, the Hamilto- 
nian can be diagonalized to the form — ^kvtvk where 
rjk are the Bogoliubov quasiparticles and Ck is the exci- 
tation energy or gap given by^li^ 

e/c = (^^ + 1 + Ji + 2hcosk - 2hJs.cos2k - 2JzC0skY/'^{A) 

with Jx set equal to unity. 

It can be easily shown that the gap of the spectrum 
vanishes at /i = J3 -I- 1 and also at /i = J3 — 1 with 
ordering (or mode-softening) wave vectors tt and re- 
spectively. These two lines correspond to quantum phase 
transitions from a ferromagnetically ordered phase to a 
quantum paramagnetic phase with the associated expo- 
nents being the same as the transverse Ising model'^S. 
The wave vector at which the minima of tk (Eq- 4) oc- 
curs, gets shifted from fc = to fc = tt wave vector when 
one crosses the line h — J3. Moreover, there is an ad- 
ditional phase transition at h = ^Jj,- This transition 
belongs to the universality class of the anisotropic tran- 
sition observed in the transverse XY-model dual to the 
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Hamiltonian (1)22, and the phase boundary is flanked by 
the incommensurate phases on either side with ordering 
wave vector given by 



cosfc — 



(5) 



This incommensurate wave vector picks up a value ko 
such that cos/cq = 1/2 J3 at the phase boundary. Obvi- 
ously, for J3 < 0.5, the anisotropic phase transition can 
not occur. The equilibrium phase diagram of the model 
is shown in Fig. 1. 



h 



Comme isurate 




FIG. 1: Equilibrium phase diagram of the three spin inter- 
acting Ising model. Solid lines show phase boundaries and 
dotted line marks the boundary between the incommensurate 
and the commensurate phase. 



III. TRANSVERSE QUENCHING AND 
RESULTS 

The dynamics of the three spin interacting TIM is 
found to be very interesting due to the fact that the 
system crosses various quantum critical lines during the 
process of dynamics. As mentioned already, the sys- 
tem deviates from the adiabatic evolution in the neigh- 
borhood of a quantum critical point where nonadia- 
baticity dominates due to the divergence of relaxation 
time. We shall introduce the time dependence in the 
Hamiltonian through the transverse field which is linearly 
quenched from —00 to -|-c» at a steady finite rate given 
by hit) ^ t/r, where the quenching time t determines 
the rate of quenching^i^i^. At time t = —00 the trans- 
verse field h — —00 and hence all the spins are pointing 
in the negative 2— direction. By virtue of the duality 
transformation, the transverse quenching of the 3 spin 
Hamiltonian corresponds to the anisotropic quenching of 
the transverse XY model where the interaction term of 
the later Hamiltonian is adiabatically changed from —00 
to o<M. 

Let us recall the Hamiltonian given in Eq. 3 with 
a time dependent transverse field h(t). This Hamil- 
tonian can be split into a sum of independent terms. 



= Efe>o-^fc(^)] where each Hk{t) operates on a 
four dimensional Hilbert space spanned by the basis vec- 
tors |0), \k,~k), |fc)and| — fc). The vacuum state where 
no c-particle is present, is denoted by |0) which corre- 
sponds to a spin configuration with all spins pointing in 
the — z-direction. The form of the Hamiltonian readily 
suggests that the parity (even or odd) of total number of 
fermions (given by nk = c\.Ck + cL^c_/c) for each mode 
is conserved. Therefore, to study the quenching dynam- 
ics, it is convenient to project the Hamiltonian Hk(t) in 
the subspace spanned by |0) and \k,—k) . The projected 
Hamiltonian has a form 



ft.(t) -I- cos fc — J3 cos 2fc i(sinfc- 
— i(sin fc — J3 sin 2fc) —{h{t) 



J3 sin 2fc) 

- cos k — cos 2fc) 



In the reduced Hilbert space, any general state can be 
represented as a superposition of |0) and |fc,— fc) with 
time dependent amplitudes Uk{t) and Vk{t) such that 
ipk{t) = Uk{t)\0) + Vk{t)\k,—k). The time evolution of 
the state is given by the Schroedinger equation 



(6) 



We shall here use the initial conditions Ufc(— 00) = 1 
and Vk{—oo) = which in the spin language corre- 
sponds to the state with all spins down. The off diag- 
onal term A = sin fc — J3 sin 2fc represents the interac- 
tion between the two time dependent levels with energies 
Ei^2 = ^[h{t) + cos(fc) — J3COs2fc]. The zeroes of the 
off-diagonal term yield the mode softening wave vectors 
fc = 0,7r and cos^^ 1/(2 J3) (provided J3 > 0.5) at which 
the system becomes quantum critical for appropriate pa- 
rameter values. At these parameter values and wave 
vectors, the system undergoes a nonadiabatic transition 
from its instantaneous ground state. A measure of nona- 
diabaticity can be obtained by comparing the two level 
problem to the corresponding Landau-Zener transition 
equations^ii. For a completely adiabatic transition, we 
expect the final state to be described by the probability 
amplitudes Ufc(-l-oo) = and Ufc(-l-oo) — 1, i.e., the com- 
plete spin-flip from down to up occurs. The nonadiabatic 
transition probability pk is directly given by |ufe(-|-oo)p 
where the probability amplitudes Uk (t) and Vk (t) are nor- 
malized at each instant of time. Equivalently, pk also 
measures the probability that the system remains in its 
initial state |0) even at the final time. Using the results 
of Landau-Zener transitions^^'^^, pk is found to be 



Pk = |wfe(+oo)p = exp(— 27r7) where 7 



Therefore, in this model 



Pk = exp[— 7rT(sin k — J3 sin 2fc)^] . 



•(7) 



(8) 



The variation of pk as a function of fc for different values 
of quenching time r is shown in Fig. 2. It is to be noted 
that for J3 < 0.5, there are peaks at — 7r,0 and tt in the 
whole range of wave vectors from — tt to tt whereas for 
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FIG. 2: Non-adiabatic transition probability pk for the three 
spin interacting Hamiltonian with J3 = 0.1 in Fig. 2(a) and 
J3 = 1 in Fig. 2(b) for various r. It should be noted that for 
J3 = 1, the system undergoes a non-adiabatic transition at 
an incommensurate wave vector k = tt/3 and therefore, there 
is an additional peak at this wave vector . For large t, pk is 
nonzero only for wave vectors very close to the critical modes. 
On the other hand, for small values of r, levels cross quickly 
resulting to a non-zero value of Pk for all values of k. 



J3 > 0.5 there are additional peaks at the incommensu- 
rate values ± cos^^(l/2 J3). 

As mentioned already, the degree of nonadiabaticity 
can be quantified through the density of kinks n gener- 
ated at t = +00 which is obtained by integrating the 
probability pk over the entire range of wave vector. 

n = ^Pfc = — y dk Pk. (9) 

A close inspection of Eq. 8 (see also Fig. 2) shows that 
for sufficiently slow quenching (i.e., large r), only modes 
close to the critical modes are excited. One can therefore, 
to the lowest order in k, replace sinfc by k in the expo- 
nential of Eq. 8 and arrive at an approximate analytical 
expression for density of kinks in the large r limit, given 
by: 



1 



n ■■ 



27r(l-2J3)V^ 27r(l + 2J3)V^ 



(for J3 < 0.5) 
(10a) 



The density of kinks monotonically increases with in- 
creasing J3 provided J3 < 0.5 because of the decrease 
in the off-diagonal term making the probability of exci- 
tations higher. On the other hand, for J3 > 0.5, the 
off-diagonal term monotonically increases with with in- 
creasing J3 resulting to an overall decrease in the density 
of kinks, see figure 3. These results can also be seen from 
the approximate analytical expression of the kink density 
given in Eq. (10 a and b) for both the cases. 

We shall now focus our attention to the case J3 = 0.5. 
In the process of the transverse quenching, the system 
crosses the multicritical point at /i = —0.5, J3 = 0.5 as 
shown in the Fig. 1 and a special power-law behavior 
of the kink density is observed at these parameter val- 
ues. The transition probability pk maximizes at fc = 
as shown above. The argument of the exponential in pk 
is expanded about k=0 at J3 = 0.5, leading to a form 
Pk = exp[— 7rTfc®/4]. The contribution to the the kink- 
density scales as l/r^/^ which can be obtained by sim- 
ply integrating this pk, see figure 3. This relatively slow 
decay of density is a special characterisitc of quenching 
through a multicritical point. A similar behavior is also 
seen in the anisotropic quenching of the transverse XY 
modeiiS. 
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FIG. 3: The variation of kink density with r for different 
J3 (< 0.5) is shown in Fig. 3a. Kink density increases with 
J3. On the other hand for J3 > 0.5, this variation decreases 
with J3 as shown in Fig. 3b. The thick line has slope —0.5 
and the slope of the dotted line is —1/6, a behavior observed 
at J3 = 0.5 



"=2.(2J3-1)V^+2.(1 + 2J3)V^ (forJ3>0.5). 

(105) 

In Eq. 10(a), the first term corresponds to the contri- 
bution from modes close to fc = whereas the second 
term is due to the peaks at k = n, —n. For the case 
J3 > 0.5, the contribution from peaks at fc = 0,7r and 
— TT happens to be the same as Eq. 10(a) whereas the 
contribution ni from the modes close to k = ko is also 
equal to 10(a) in the following way: 

1 r . 

m — — exp^ 



[— 7rT{{cos A:q — 2J3 cos 2fco)(^'~^'o)}^] 



1 



1 



1 



27r/?'2J3 + l 2J3-1^ 



(11) 



One can also study the effect of the anisotropic quench- 
ing which involves quenching of the nearest neighbor 
Ising interaction term J^(i)(~ f/r) from —00 to +00 
instead of the transverse field with the three spin inter- 
action term set to unity. At i — > —00, the ground state of 
the system is antiferromagnetic along x. The probability 
of the non-adiabatic transition is similarly given by 

Pfe = exp[-7rr{(/i -I- 1) sinA:}^]. (12) 

It is interesting to note that for h = —1, pk is unity for 
all values of k. The density of kinks for the anisotropic 
quenching is given as 



n = — dkpk — exp((— 7rr(/i + 1) sin/c)^ 
27r J 



(13) 
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with an approximate analytical form given as 
1 



(14) 



which shows that the density of kinks decreases mono- 
tonically with h. This can be attributed to an increase 
in the off-diagonal term of the Hamiltonian. 



IV. CONNECTION TO THE TRANSVERSE 
QUENCHING OF THE MODELS WITH NEXT 
NEAREST NEIGHBOR INTERACTIONS 

We shall now use the results of the previous section 
to study the transverse quenching of a quantum Ising 
model with uniform ferromagnetic nearest neighbor in- 
teraction and also an additional NNN interaction which 
is either antiferromagnetic or ferromagnetic. The model 
with NNN antiferromagnetic interaction has regular frus- 
trations and is popularly known as Axial Next Nearest 
Neighbor Ising (ANNNI) model^ in a transverse field. 
We shall show below that within a mean field approxi- 
mation, the three spin model has a close resemblance to 
the one-dimensional NNN interacting TIMs. 

The Hamiltonian of the transverse ANNNI model is 
given by 



1 ^ 



(15) 



where Ji, J2 > 0. Henceforth, without loss of any gen- 
erality, we shall set Ji = 1. At /i = 0, the ground state 
is ferromagnetically ordered for J2 < 0.5, whereas the 
system shows an "anti-phase" ordering (where two up 
spins are followed by two down spins) for J2 > 0.5. The 
two phases meet at an infinitely degenerate multi-critical 
point J2 = 0.5 and h = Q. The quantum fluctuations 
introduced by the transverse field h competes with the 
ferromagnetic (or the antiphase) order and eventually 
the system undergoes a quantum phase transition to a 
paramagnetic phase at a critical value of the transverse 
field given by he which is a function of the NNN interac- 
tion J2. One-dimensional quantum ANNNI model shows 
a rich phase diagram which is not fully understood till 

When mapped to the corresponding fermionic Hamil- 
tonian via a JW transformation, the NNN interaction 
leads to a four-fermion term in the fermionic version of 
the Hamiltonian. In the limit J2 — > 0, this term van- 
ishes so that the model is exactly solvable in terms of 
non-interacting fermions. For non-zero J2, the fermionic 
Hamiltonian is written as 

H^^\{[Y. ^^(24^. - 1) + (cl - cO(cl+i + c.+i) 

i 

-J2{c\ ~ q)(1 - 2c|+iQ+i)(4+2 + C..+2)]}. (16) 

The occurrence of the four-fermion term renders the 
model analytically intractable though an approximate 



analytical solution is possible at least in the limit of small 
J2. Deep in the paramagnetic phase all the spins are 
oriented in the direction of the transverse field so that 
< erf >= 1 or in the fermionic language 1 — 2c|^j^Ci+i = 
— 1. We shall approximate < erf >= 1 for all pos- 
itive values of h including h ^ Q. This approxima- 
tion, though crude, transforms the four fermion term 
{c\ - Cj)(l - 2c|_,_iCi+i)(c-_,_2 -I- Ci+2) to a quadratic form. 
The Hamiltonian becomes exactly solvable but the rich 
phase diagram of the model is not captured in this 



approximation 



35,36 



Within this approximation, we shall 



explore the role of small NNN antiferromagnetic interac- 
tion on the density of kinks produced during the trans- 
verse quenching. As described below, this approximation 
at least shows a decrease of critical field he with J2 for 
J2 < 0.5. 

The mean field Hamiltonian in the momentum space 

is 

H = — [^(/i -I- cosA: -I- J2Cos2A:)(cj,Cfe -|- ctj,c„fc) 



k>0 



-|-i(sinA: + J2sin2fc) (c|.c^j, -f CkC-k)] (17) 



Comparing Eq. (17) with Eq. (3), one finds that the 
transverse ANNNI chain Hamiltonian in the mean field 
approximation is identical to the three spin Hamiltonian 
if the antiferromagnetic interaction J2 of the former is re- 
placed by the negative of the three spin interaction term 
(J3) in the latter. Using the results of the previous sec- 
tion, the phase diagram of the mean field ANNNI model 
can be found out for h > 0(see Fig. 4). The phase bound- 
ary between the ferromagnetic phase and the paramag- 
netic is given by h — 1 — J2 with an ordering wave vector 
TT (this corresponds to the Ising transition at h = J3 + 1 
of Fig. 1). For J3 > 0.5, i.e., the transition between the 
antiphase and the paramagnetic phase, is given by the 
corresponding anisotropic transition of three spin model 
with the phase boundary given by the equation h = J2 
and the ordering wave vector has an incommensurate 
value as given in Eq. 5. 

The approximation l~2c\ci ~ — 1 is valid for positive h 
only, we choose a quenching scheme where the transverse 
field has a functional dependence h{t) ~t/T with t 
going from —00 to so that h(t) remains positive for the 
entire quenching period and vanishes at the end of the 
quenching. Therefore the system does not cross the Ising 
critical line —h= J2 + 1. 

In the final state at t = 0, all the spins are expected to 
orient in the x-direction with a ferromagnetic order. The 
density of oppositely oriented spins at f is related to 
J2 as 



1 



27r(l - 2J2)V7 



(18) 



which shows that the density of kinks increases monoton- 
ically with J2. 

It should be noted that if we follow a quenching scheme 
in which the transverse field is changed from —00 to zero. 
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Quenching direction 




Fen'omagnet Antiphase 



0.5 



FIG. 4: Mean field phase diagram of the ANNNI model in 
the h — J2 plane. We study the quenching dynamics across 
the phase boundary close to J2 —>■ 0. 



with +1 



we must approximate the term 1 — 2c\^^CiJ^i 
rather than —1 for above calculations to be viable. In 
the process of dynamics, the system crosses the quantum 
critical line —h~ J2 + 1 with the modes close to fc = 
getting excited. This approach as well leads to identical 
result for kink density (as given in Eq. 18). Therefore, 
the presence of a small antiferromagnetic NNN interac- 
tion adds to the kink-production in comparison to the 
ferromagnetic transverse Ising model (J2 = 0) with the 
same quenching scheme. 

One can also study, in the similar spirit, a model with 
a small ferromagnetic NNN interaction JpM- We use 
the same mean field approximation for /i > so that this 
model is identical to the three spin model with J3 = Jfm- 
A similar calculation leads once again to a 1/ \/t fall of 
the density of kinks given by 



1 



(19) 



This is expected because the NNN ferromagnetic inter- 
action enhances the strength of the ferromagnetic order- 
ing and hence the probability of excitations or density of 
kinks is lowered. 



CONCLUSIONS 



from model to model. The first of the variants includes a 
three spin interaction with strength J3. Here, the phase 
diagram indicates the existence of an anisotropic phase 
transition at an incommensurate value of wave vector in 
addition to the normal Ising transition for J3 > 0.5. In- 
terestingly, we observe that the density of kinks increases 
monotonically with J3 for J3 < 0.5 whereas decreases for 
J3 > 0.5. On the other hand, at J3 = 0.5, the con- 
tribution to the kink density scales as r~^/^ due to the 
existence of a multicritical point at J3 = 0.5. The other 
set of Hamiltonians include a ferromagnetic or an anti- 
ferromagnetic next nearest neighbor interactions. The 
presence of the four fermion term makes such a Hamil- 
tonian analytically intractable. We have used a mean 
field approximation to reduce the four fermion term in 
the fermionized representation to a quadratic term. The 
quenching scheme is chosen carefully so that the regions 
where the approximation is not valid are avoided in the 
process of dynamics. Using the similarity between the 
fermionized next nearest neighbor interacting Hamilto- 
nians under the mean field approximation, and the three 
spin interacting model, the density of kink in the final 
state is estimated. It is observed that the ferromagnetic 
next nearest neighbor interactions reduces the density 
of kinks produced as opposed to the case of antiferro- 
magnetic next nearest neighbor interaction because such 
a ferromagnetic interaction enhances the ferro-ordering 
discouraging the production of kinks. On the other hand, 
frustration leads an enhanced non-adiabatic transitions. 
We should mention in conclusion that it is in principle 
possible to construct a better mean field theory for the 
ANNNI model^, however, no qualitative change in the 
dynamical behaviour in the region J2 — > is expected. 

We conclude with the comment that the models stud- 
ied in the present work are integrable (at least in the 
mean field limit) which leads to a 1/v^ scaling behav- 
ior of the defect density. However, in a random or 
a non-integrable system such a behavior need not be 
expected^*^. The quenching and annealing dynamics of 
several non-integrable systems along with the depen- 
dence of the defect density on the integrability of the 
model are yet to be completely understood. We have also 
observed a much slower decay of the form \/t^/^ when 
quenched through the multicritical point of the three spin 
model as in the anisotropic quenching of the transverse 
XY chainlS. 



In this paper, we have studied the effect of various 
additional interactions on the dynamics of the transverse 
Ising model when swept across the quantum critical lines. 
The defect density scale with the timescale r as r"-'^/^, 
like in transverse Ising case, with a prefactor which varies 
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